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I use the recently proposed framework of ontological models [Harrigan et ai, arXiv:0709.1149v2] 
to obtain economical models for results of tomographically complete sets of measurements on finite- 
dimensional quantum systems. I describe a procedure that simplifies the models by decreasing the 
number of necessary ontic states, and present an explicit model with just 33 ontic states for a qutrit. 
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I. INTRODUCTION 

The inherent complexity in the description of quan- 
tum states and operations must be studied if we are to 
understand the foundations of the theory. This study 
can shed Ught on the use of quantum systems for infor- 
mation processing, and clarify issues of simulabiUty and 
compatibiUty with other theories, such as relativity. 

In this context, there have been some recent proposals 
to describe quantum theory using general probabiUstic 
theories. The operational "r-p" framework 0,0, al- 
lows us to study general probabilistic theories, of which 
quantum mechanics is a particular instance. With the 
interest in quantum information and in simulating com- 
plex quantum phenomena, these recent approaches have 
yielded interesting results bearing on the efficiency with 
which we can simulate certain classes of quantum op- 
erations 0, Q. They have also contributed to the un- 
derstanding of negativity and contextuality in quasi- 
probability distributions ^TJ. 

Harrigan et al. have recently proposed the onto- 
logical model formalism to describe any compilation of 
experimental outcomes. This framework is based on re- 
alism, and involves only positive probability distributions 
and indicator functions. I review their formalism in sec- 
tion [TTl They presented three types of general models 
that describe any set of experimental outcomes. In their 
paper, however, the models are applied only to a few 
quantum examples that involved discrete sets of prepa- 
rations (states) and measurements. A similar approach 
was recently used to obtain asymptotically efficient mod- 
els for general quantum measurements [9] . 

In section IIII Al I apply Harrigan et al.^s framework 
to obtain a model for data arising from a tomograph- 
ically complete set of measurements on any state of a 
d-dimensional quantum system. I describe how one may 
decrease the complexity of the model (as measured by the 
number of necessary ontic states), and obtain a model for 
a qutrit {d = 3) with just 33 ontic states. Finding such 
economical models is important for problems such as the 
study of hidden-variable theories for quantum mechan- 
ics, and the search for efficient descriptions of general 
experimental data (see Q for further references). 



II. THE ONTOLOGICAL MODEL FORMALISM 

In this section I review the ontological model formal- 
ism, introduced by Harrigan, Rudolph and Aaronson in 
Q . The aim of the formalism is to obtain a model that 
describes experimental data compiled in a data table. 
The model uses positive, normalized probability distri- 
butions over a certain number of variables, and indicator 
functions to describe measurements. 

A data table is a compilation of the (probabilistic) 
results of a series of experiments. An experiment con- 
sists of a set of s preparation procedures V^^\ where 
i = l,2,...,s, together with a set of m different mea- 
surement procedures where i = 1,2, m. We can 
choose the number of outcomes for each measurement to 
be d (i.e. the same for all measurements), if necessary by 
padding with null outcomes. 

For concreteness, let us see how to describe a data table 
arising from projective measurements on a qubit. Each 
preparation represents a quantum state. The measure- 
ments we will consider are the set {A^^^^ = X,Ai'--^^ = 
Y,M'^"^^ = Z} (the Pauli matrices). This set of mea- 
surements is tomographically complete, meaning that 
the probabilities associated with these three operators 
are sufficient to completely characterize the preparations, 
and hence sufficient also to calculate the probabilities as- 
sociated with outcomes of any other measurement. 

After many measurements of observables the 
probability (relative frequency) of each outcome can be 
tabulated in a data table. Below we represent three data 
tables, corresponding to each of the three measurements 
For convenience, the three data tables are repre- 
sented on top of each other, forming a single, composite 
data table: 
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In the table each column corresponds to a preparation 
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T-'*'-'^. The probabilities should match those computed 
from quantum theory via the trace rule. Note that since 
there is a continuum of pure quantum states, our table 
has an infinite number of columns. I have kept the nota- 
tion simple by pretending that the index parameterizing 
the preparations is discrete. It is clear that the table rep- 
resents the two probabilistic outcomes of measurements 
(rows), which will of course depend on the preparation 
used (columns). 



ourselves with a simple ontological compression scheme 
for data arising from a particular tomographically com- 
plete set of measurements on c?-dimensional systems. 
We will start with Harrigan et aZ.'s general model with 
— ontic states (which is valid even for s = oo), and 
obtain economical OF's by carefully excluding some of 
the original ontic states. Let us start by analyzing the 
simplest such case, which consists of measurements on a 
2-level system, i.e. a qubit. 



A. Ontological factorizations of a data table 

We start with a composite data table D correspond- 
ing to TO d— outcome measurements. The probabilis- 
tic outcomes are compiled in m c? x s matrices D^^\ 
X ~ 1, 2, ...,TO. An ontological factorization of data ta- 
ble D Q consists of a factorization of each D^^^ into the 
product of a d X O measurement matrix A^^^^ and an 
il X s preparation matrix V: 

=>l(^)7', Va; = 1,2,...TO. (2) 

The matrices Al^^^ and V must be non-negative and 
column-stochastic, that is, < Vjk < 1, and Vjk = 1 
for all J (and similarly for MS^^). 

Ontological factorizations are not unique and represent 
models, having the following interpretation. The fcth col- 
umn of V represents a classical probability distribution 
over f2 ontic states, which corresponds to the preparation 

procedure V'^^'' of D. As for tW^'^^ ^^ij '^'^ proba- 
bility of obtaining measurement outcome z, if the ontic 
state is j, when performing measurement x. Note that 
the complexity of the model (measured by the size of the 
matrices involved) increases with the number il of ontic 
states underlying it. 

If all matrices M. have entries which are either or 1 
the ontological factorization is called deterministic^ oth- 
erwise it is indeterministic. In Q Harrigan et al. ob- 
tained a few general constructions of ontological factor- 
izations for data tables. Let us call the number of mea- 
surements TO, the number of preparations s, and the num- 
ber of outcomes of each measurement d. Then in [8] three 
OF constructions were presented: an indeterministic OF 
with f2 = s ontic states; a deterministic OF with = d™ 
ontic states; and a more parsimonious (on to) determin- 
istic OF with 17 — s{dm — l) ontic states. These OF's are 
general, that is, they represent valid models for arbitrary 
data tables. 

OF's can be said to be economical when they have 
fewer ontic states than the general models above. Eco- 
nomical models decrease the number of ontic states by ex- 
ploring regularities in the data tables. In the case treated 
in the next sections, these regularities are linear inequal- 
ities satisfied by the probabilities arising from quantum 
measurements. 

The process of obtaining economical OF's was dubbed 
ontological compression. In what follows we will concern 



B. Example: one qubit 

In this section I illustrate the ideas reviewed above by 
presenting an ontological factorization (OF) for the data 
table (dl), which refers to the set of all preparations for 
a single qubit, and measurements of the Pauli operators 
X,Y, and Z. I chose to present this data table and its 
OF's here as an example to be followed for larger systems 
in the next section. 

We could have chosen to study more general data ta- 
bles for a qubit, including for example the continuum of 
possible projective measurements, as opposed to just the 
three Pauli operators. It is well known, however, that 
the probabilities associated with just three independent 
two-outcome observables are sufficient to calculate the 
probabilities associated with any of any other observ- 
able. Moreover, the eigenvectors of the Pauli operators 
form mutually unbiased bases, which are optimal for to- 
mography. 

Harrigan et aZ.'s deterministic OF for this data table 
has = 8 ontic states. The motivation for this particular 
OF is as follows. Experimenting with different prepara- 
tions, one finds that there exist states yielding definite, 
repeatable results for each of the three measurements. 
These are, of course, the eigenstates of the X, Y, Z ob- 
servables. The most general model for three two-outcome 
experiments would allow for simultaneous definite out- 
comes for these three observables, which of course is not 
allowed by quantum mechanics. Nevertheless, one can 
understand quantum mechanics as a constraint on these 
more general probabilistic models, an idea that has been 
around for a long time P, d, |3|, \4, llfll ■ 

With this in mind, the set of ontic states I propose 
consists of all possible definite states for these three ob- 
servables. Since each observable has two outcomes, there 
are 2"^ = 8 such ontic states. We can represent these on- 
tic states by the columns of the 8x8 identity matrix; it 
is clear that arbitrary probability n-tuples can be writ- 
ten as convex combinations of these 8 8-tuples. We need 
to assign a meaning to the position of the 1 in each on- 
tic state, before we can present the matrices representing 
measurements. Thus, let j be the position of the 1 in 
an ontic state (j — 1, 2, 8). In binary, j — 1 is written 
as a three-bit string 636261, with bit hi indicating which 
of the two outcomes (labeled or 1) should occur for a 
measurement of basis i. For example, j = 5 corresponds 
to the ontic n-tuple P5 = (0,0,0,0,1,0,0,0)^. The bi- 
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nary representation of j — 1 = 4 is 100, indicating that 
this ontic state would yield definite outcomes of 0,0 and 
1 for measurements of first, second and third bases, re- 
spectively. 

The matrices representing these measurements were 
chosen according to the convention I used in numbering 
the ontic states. The three measurements are represented 
as follows: 




10 10 10 10 

10 10 10 1^ 

1 1 1 1 o' 

1 1 1 1^ 

1 1 1 1 
1 1 1 1 



(3) 
(4) 
(5) 



As an illustration, suppose we want to measure Z on 
a state (preparation) represented by the following one- 
column matrix V: 



(6) 



where Vi is the probability of preparation in ontic state 
i. The data table Z?*^^^ with the probabilities associated 
with outcomes (po) and 1 (pi) of Z is given by the 
factorization 



Pi 



(7) 



Any data table D^^> corresponding to three two- 
outcome measurements can be given an OF using these 
eight ontic states - this is an example of one of the general 
constructions proposed in Q . For m d-outcome measure- 
ments and s preparations, the number of necessary ontic 
states is at most 51 = d™, but this number can be smaller, 
depending on the data table. We are interested in finding 
economical OF's for data arising from quantum measure- 
ments. 

Since the probabilities associated with the two out- 
comes of each observable must add up to 1, we can rep- 
resent experimental outcomes by listing just one indepen- 
dent probability per observable, say the first one. For ex- 
ample, the first column in table ([1]) can be represented by 
the 3-tuple p = (0.93,0.73,0.59). In this 3-dimensional 
space of probabilities, general preparations are the con- 
vex hull of the eight ontic state vertices, which form a 
0/1 cube in 3 dimensions which we will call polyhedron 
Pi. Quantum states are a proper subset of the cube, the 



1/2 and with center p 
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ball of radius r 

Fig. 

To perform ontological compression, we need to find a 
smaller set of ontic states whose convex combinations still 
account for the data table. The simple ontological com- 
pression scheme I will consider consists of keeping only 
a proper subset of the original model's ontic states, and 
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Pi 



P2 (b) 



FIG. 1: (Color online) (a) The cube Pi represents general 
probability 3-tuples, which include the ball representing quan- 
tum states of a qubit. (b) If we eliminate the cube vertex 
(1, 1, 1), the convex hull of the remaining 7 vertices is polyhe- 
dron P2, which no longer contains all quantum states. All 8 
vertices are required to account for all quantum states in this 
way. 



adapting the measurement matrices accordingly. For ex- 
ample, if we eliminate the eighth ontic state (cube vertex 
p — (1, 1, 1)), the allowed preparations now can be con- 
vex combinations of only the 7 remaining ontic states. 
Geometrically, this is the convex hull of 7 vertices, the 
polyhedron P2 represented in Fig. [T]-(b). This deletion 
of one ontic state must be accompanied by a correspond- 
ing change in the measurement matrices ^A^^\M.^'^^ and 
M.^'^'^ - the column associated with the deleted ontic state, 
the eighth, must also be deleted. 

By comparing the new polyhedron P2 with the ball of 
quantum states, we see that there exist quantum states 
p ^ P2- Hence, the proposed ontological compression 
scheme fails. Symmetry indicates the compression will 
also fail in case of removal of any other ontic state from 
our model. In this case we need the full set of eight ontic 
states if we are to represent all quantum states as convex 
combinations of a subset of the original set of ontic states 
(the cube vertices). 

In the next section I will show that this failure of the 
simple ontological compression scheme I propose happens 
only for d = 2. For d > 3 the scheme outlined above 
can greatly reduce the number of ontic states necessary 
to describe data tables arising from a tomographically 
complete set of measurements. 



III. ECONOMICAL ONTOLOGICAL 
FACTORIZATIONS 

In this section I show how to obtain ontological factor- 
izations for data tables associated with all preparations 
of a d-dimensional quantum system, together with mea- 
surements onto a particular tomographically complete set 
of observables. I will start with the general, extravagant 
OF with n — d"^ ontic states of Harrigan et al. [1] , and 
show how the number of ontic states can be reduced, 
resulting in economical OF's. 
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The measurements we will consider are projectors onto 
a complete set of mutually unbiased bases (MUB) for a 
d-dimensional system. Such a set contains d + 1 bases, 
each providing d — 1 independent measurement proba- 
bilities. Such complete sets of MUB are known to ex- 
ist for any Hilbert space dimension d which is power-of- 
prime. In the case of prime-dimensional systems, there is 
a canonical construction of a complete set of MUB, due 
to Ivanovic [ll| - in the Appendix I list the four bases 
for d = 3. We will use this canonical set, and restrict 
ourselves to discussing the case of prime d. 

Our initial model is Harrigan et aL's deterministic OF 
Q, in which each measurement can yield any definite 
result. This OF must have as many ontic states as there 
are definite outcomes for the d + 1 MUB measurements. 
Hence, the number of ontic states for this OF is d'^'^^. We 
will denote the polytope formed by taking the convex hull 
of these ontic states Pi. The measurement matrices can 
be defined in a similar way as we did for a qubit in the 
last section, resulting in a deterministic OF. 

In order to find OF's with fewer ontic states, we resort 
to some ideas from convex geometry. The ontic 
states can be chosen so that the n-tuple representing the 
jth ontic state has O's in all positions, except the jth, 
which has 1. The set of allowed preparations are convex 
combinations of those ontic states, i.e. polytope Pi. If 
we remove ontic states from that initial set, the resulting 
OF will only account for preparations which lie on the 
convex hull of the remaining ontic states. This convex 
hull is a polytope, which we will call P2. Our simple 
ontological compression scheme will involve finding sub- 
sets of the original set of ontic states which include all 
quantum states in their convex hull. 

The identification of OF's with fewer ontic states is a 
way of approximating the convex set of quantum states 
by simpler convex bodies, the polytopes. The facets 
of the enclosing polytopes represent linear inequalities 
obeyed by any quantum state. These inequalities are 
the regularities we identify in the quantum data tables, 
enabling us to do the ontological compression. In this 
representation using probabilities associated with MUB, 
all pure quantum states lie on the surface of a high- 
dimensional ellipsoid 

Different approaches can be taken to perform this sim- 
ple kind of ontological compression. In what follows I 
describe a general procedure that works for any d > 3, 
and exemplify the approach by tackling the qutrit case 
(d = 3). 



A. A general compression scheme 

Here I will point out one way of obtaining economical 
ontological factorizations (OF's) for d-dimensional sys- 
tems, for prime d. First, we need to review how we can 
represent states using probabilities associated with mu- 
tually unbiased bases (MUB). 

Not all components of the probability n-tuple p are 



equivalent: sets of d — 1 components belong together, as 
they are independent probabilities associated with dif- 
ferent outcomes of the same d-outcome observable. Let 
greek indices k, ^ etc denote the observable associated 
with a given set of probabilities. The representation of 
states in probability space would be with a (d^ — l)-tuple 

of probabilities p'j^^' 

p = (p(^p(^...,p«,,...,pr\...,p[j^i;), (8) 

where the super-indices indicate which of the MUB k the 
probabilities are associated with. Our initial OF consists 
of measurement matrices M that are trivial generaliza- 
tions of the qubit case, and of the d'^+^ ontic states rep- 
resented by probability vectors as above, where for each 
of the MUB K either all p^"'* = (indicating that the 
omitted p^''' — 1); or there is exactly one of the pj"' = 1 

for each k, with the remaining p['^'^ — 0. The convex hull 
of this set of d'^^^ ontic states is our initial polytope Pi. 
This choice of initial ontic states results in a polytope in 

(k) 

which the probabilities p] obey the MUB constraints 

= 1 (for aU k). 

With a view to obtaining OF's which are as economical 
as possible, we can start with polytope Pi and remove 
one randomly chosen vertex (a process that we can later 
iterate). The facets of the smaller resulting polytope P2 
can be found with a convex hull algorithm such as the 
one presented in [l^. To guarantee that all quantum 
states are contained in P2, it is enough to consider each 
facet in turn, proving that no quantum state violates it. 
This can be done as follows. 

For simplicity of notation, let us use a single sub-index 
i to index the {d^ — 1) independent probabilities aris- 
ing from MUB measurements that fully describe a d- 
dimensional quantum state. Each probability pi is ob- 
tained from the associated MUB projector Pi using the 
trace rule (eq. lA.ip . 

Each facet is represented by a linear inequality on the 
MUB probabilities, and as such can be written as the ex- 
pectation value of an operator built out of the associated 
projectors. The association is as follows: 

^ C^P^ > / ^ ^ /' 

where / is a constant offset. Because all coefficients Ci 
are real, the operator C = y^.- c,:P,: is Hermitian. In this 
case, the minimax principle ([3]) §90) guarantees that 
the minimum quantum expectation value of C will be 
given by its smallest eigenvalue Xmin- If ^mm > / then 
no quantum state falls outside the facet. 

B. Example: d = 3 

For d = 3, I started with the 81-vertex polytope Pi 
and tried to remove vertices randomly, checking whether 



5 



1 



5 



10 



15 



20 



25 



30 




13 5 7 



FIG. 2: The rows represent the 33 vertices of an 8-dimensional 
polytope that includes aU qutrit states, where black and while 
code for 1 and 0, respectively. 



all quantum states satisfy the inequalities representing 
the facets of the resulting polytope. This was done with 
the convex hull algorithm of [l^ and simple eigenvalue 
calculations as described in the last section. 

By iterating this procedure many times (backtrack- 
ing when necessary), I obtained economical, determin- 
istic OF's for a qutrit with as few as 33 ontic states. In 
Fig. [2] I list the vertices of one such OF. The probabili- 
ties are associated with the MUB projectors listed in the 
Appendix. 

This 33-vertex polytope can be described also in terms 
of its facets. These are represented by 51 linear inequal- 
ities on the probabilities. Of these, 12 are the trivial 
inequalities X^iPi"'' ^ 1 (4 inequalities corresponding to 
K = 1,2,3,4) and > (8 inequalities corresponding to 
K — 1,2,3,4 and j = 1, 2 ). The remaining 39 inequalities 
are listed as the rows of Table HI Each inequality is of the 
form Y^^=i '^iPi — /■ The last column of Table |T] presents 
the smallest eigenvalue Xmin of operator C = ciPi, where 
Pi is the projector associated with probability pi. By the 
results of the last section, we see that no quantum state 



TABLE I: Polytope facets 
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C2 


C3 


C4 


C5 


C6 


C7 


C8 
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-2 


-1 


-2 


-1 


1 


-1 


-1 


-2 


-5 


-4.3028 
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-1 


2 
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-2 


-1 


-4 


-3.8608 
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-1 


-1 
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-2 
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-3 


-2 
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1 


-2 


-3 


-5 


-4.8455 


-1 


-2 


-1 
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-2 


-1 


-2 


-1 


-6 
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-4 
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1 


-1 
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-2.0000 
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2 


1 


-1 


2 


2 


-1 


1 
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2 
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-3 


-2 


-3 


-1 


-2 


-3 


-7 


-6.8455 


2 


1 


-2 


-1 


-1 


1 


-1 


-2 


-3 


-2.8608 



falls outside the polytope, and also that some of the facet 
inequalities are saturated by quantum states. At least 
some of the inequalities presented in Table [J are new, see 
[l5| for a different set of similar linear inequalities, de- 
rived in finding the extrenia of discrete Wigner functions 
for prime dimension d. 

I have also checked that this 33-vertex polytope is min- 
imal, in the sense that removal of any of its 33 vertices 



6 



results in a polytope that no longer includes all quan- 
tum states in its interior, i.e. its ontic states no longer 
account for all quantum states as convex combinations. 
There are many other minimal polytopes (with varying 
number of vertices) representing economical OF's for a 
qutrit. A limited random search failed to find any such 
OF's with less than 33 ontic states. 

This computational approach to finding economical 
OF's can complement other approaches in the search for 
ontological models with as few ontic states as possible. 
For example, one may start with the construction of on- 
tological models recently delineated in Q and apply the 
random pruning described here to decrease the number 
of ontic states necessary in an ontological model for to- 
mographically complete sets of quantum measurements. 



IV. CONCLUSION 

I have reviewed the ontological model formalism pro- 
posed in [s], and described a procedure that results in 
economical models for a particular set of tomographi- 
cally complete measurements on general states of a d- 
dimensional quantum system. To illustrate the proce- 
dure, I brought down the number of necessary ontic states 
from 81 to just 33 in the case of a qutrit {d = 3). In the 
process, I found some new linear inequalities for probabil- 
ities associated to mutually unbiased bases for a qutrit. 
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APPENDIX 



In this Appendix I list the four mutually unbiased 
bases (MUB) for a qutrit that we used in section IIIIBI 

The probabilities pl'^'^ of measuring density matrix p onto 
the ith projector of basis k are defined as: 



(A.l) 



Using uj = exp (27rz/3), the MUB vectors can be writ- 
ten as: 





= (0,1,0) 


(A.2) 




= (0,0,1) 


(A.3) 




= (1,0,0) 


(A.4) 
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- l/x/3(l,l,l) 


(A.5) 
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= 1/Vs{lo,1,lu^) 


(A.6) 
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v\ 


= l/\/3(w2,l,w) 


(A.7) 


(3) 

v\ 


= 1/V3(w,l,l) 


(A.8) 


(3) 


- l/\/3(l,c^,l) 


(A.9) 


(3) 


= 1/V3(l,l,c^) 


(A.IO) 


(4) 


^ 1/V3{lj\1,1) 


(A.ll) 




= 1/a/3(1,c.2,1) 


(A.12) 


4'^ 


= 1/V3(l,l,c^2) 


(A.13) 



This is the d = 3 case of the canonical construction 
of complete sets of mutually unbiased bases for prime- 
dimensional quantum systems due to Ivanovic 
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